Solar Physics 
DOI: 10.1007/' 



Segmentation of Loops from Coronal EUV 
Images 

B. Inhester^, L. Feng^'^, and T. Wiegelmann^ 

Received: 26 April 2007, accepted: July 31 2007 
doi: 10.1007/sll207-007-9027-l 

© Springer •••• 

Abstract We present a procedure which extracts bright loop features from solar 
EUV images. In terms of image intensities, these features are elongated ridge-like 
intensity maxima. To discriminate the maxima, we need information about the 
spatial derivatives of the image intensity. Commonly, the derivative estimates are 
strongly affected by image noise. We therefore use a regularized estimation of the 
derivative which is then used to interpolate a discrete vector field of ridge points 
"ridgels" which are positioned on the ridge center and have the intrinsic orientation 
of tfie local ridge direction. A scheme is proposed to connect ridgels to smooth, 
spline-represented curves which fit the observed loops. Finally, a half-automated 
user interface allows one to merge or split, eliminate or select loop fits obtained 
form the above procedure. In this paper we apply our tool to one of the first EUV 
images observed by the SECCHI instrument onboard the recently launched STEREO 
spacecraft. We compare the extracted loops with projected field lines computed from 
almost-simultaneously-taken magnetograms measured by the SOHO/MDI Doppler 
imager. The field lines were calculated using a linear force-free field model. This 
comparison allows one to verify faint and spurious loop connections produced by 
our segmentation tool and it also helps to prove the quality of the magnetic-field 
model where well-identified loop structures comply with field-line projections. We 
also discuss further potential applications of our tool such as loop oscillations and 
stereoscopy. 

Keywords: EUV images, coronal magnetic fields, image processing 
1. Introduction 

Solar EUV images ofl^er a wealth of information about the structure of the solar chro- 
mosphere, transition region, and corona. Moreover, these structures are in continuous 
motion so that the information collected by EUV images of the Sun is enormous. 
For many purposes this information must be reduced. A standard task for many 
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applications, e.g., for the comparison with projected field lines computed from a 
coronal magnetic-field model or for tie-point stereoscopic reconstruction, requires 
the extraction ofthe shape of bright loops from these images. 

Solar physics shares this task of ridge detection with many other disciplines in 
physics and also in other areas of research. A wealth of different approaches for 
the detection and segmentation of ridges has been proposed ranging from multiscale 
filtering (KoUer et al., 1995; Lindeberg, 1998) and curvelet and ridgelet transforms 
(Starck, Donoho, and Candes, 2003) to snake and watershed algorithms (Nguyen, 
Worring, and van den Boomgaard, 2000) and combining detected ridge points by 
tensor voting (Aledioni, Tang, and Lee, 2000). These general methods however al- 
ways need to be modified and optimized for specific applications. Much work in this 
field has been motivated by medical imaging {e.g., Jang and Hong, 2002; Dimas, 
Scholz, and Obermayer, 2002) and also by the application in more technical fields 
like fingerprint classification (Zhang and Yan, 2004) and the detection of roads in 
areal photography (Steger, 1998). 

For the automated segmentation of loops, a first step was made by Strous (2002, 
unpublished) who proposed a procedure to detect pixels in the vicinity of loops. This 
approach was further extended by Lee, Newman, and Gary (2006) by means of a 
connection scheme which makes use of a local solar-surface magnetic-field estimate 
to obtain a prefereable connection orientation. The procedure then leads to spline 
curves as approximations for the loop shapes in the image. The method gave quite 
promising results for artificial and also for observed trace EUV images. 

The program presented here can be considered an extension of the work by Lee, 
Newman, and Gary (2006). The improvements which we propose are to replace 
Strous' ridge-point detection scheme by a modified multiscale approach of Linde- 
berg (1998) which automatically adjusts to varying loop thicknesses and returns 
also an estimate of the reliability of the ridge point location and orientation. When 
connecting the ridge points, we would like not to use any magnetic field information 
as this prejudices a later comparison of the extracted loops with field lines computed 
from an extrapolation of the surface magnetic field. As we consider this comparison 
a validity test for the underlying field extrapolation, it would be desirable to derive 
the loop shapes independently. Our connectivity method is therefore based only on 
geometrical principles which combines the orientation of the loop at the ridge point 
with the cocircularity constraint proposed by Parent and Zucker (1989). 

The procedure is performed in three steps, each of which could be considered a 
module of its own and performs a very specific task. In the following sections we 
explain these individual steps in some detail. In the successive section we apply the 
scheme to one of the first images observed by the SECCHI instruments onboard the 
recently-launched STEREO spacecraft (Howard et al., 2007) in order to demonstrate 
the capability of our tool. Our procedure offers alternative subschemes and adaptive 
parameters to be adjusted to the contrast and noise level of the image to be dealt 
with. We discuss how the result depends on the choice of some of these parameters. 
In the final section we discuss potential applications of our tool. 

2. Method 

Our approach consists of three modular steps, each of which is described in one of the 
following subsections. The first is to find points which presumably are located on the 
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loop axis. At these positions, we also estimate the orientation of the loop for these 
estimates. Each item with this set of information is called a ridgel. The next step is 
to establish probable neighbourhood relations between them which yields chains of 
ridgels. Finally, each chain is fitted by a smoothing spline which approximates the 
loop which gave rise to the ridgels. 

2.1. Ridgel Location and Orientation 

In terms of image intensities, loop structures are elongated ridge-like intensity max- 
ima. To discriminate the maxima, we need information about the spatial derivatives 
of the image intensity. Commonly, these derivatives are strongly affected by image 
noise. In fact, numerical differentiation of data is an ill-posed problem and calls for 
proper regularization. 

We denote by i £ the integer coordinate values of the pixel centres in the 
image and by a; G the 2D continuous image coordinates with x = i at the pixel 
centres. We further assume that the observed image intensity varies sufficiently 
smoothly so that a Taylor expansion at the cell centres is a good approximation to 
the true intensity variation I{x) in the neighbourhood of i, i.e., 

I{x) ~ i{x) =c + g'^{x-i) + {x- ifH{x - i) (1) 

Pixels close to a ridge in the image intensity can then be detected on the basis of 
the local derivatives g and H (the factor 1/2 is absorbed in H). We achieve this by 
diagonalizing H, i.e., we determine the unitary matrix U with 

t/^^i/f/ = diag(/ij_, /i II ) where U = {u_\_,u\\) (2) 

where we assume that the eigenvector columns ttj^ and tt|| of U associated to the 
eigenvalues h± and /i||, respectively, are ordered so that h± < Zip. 

We have implemented two ways to estimate the Taylor coefficients. The first is a 
local fit of (1) to the image within a (2m + 1) x (2m + 1) pixel box centered around 
each pixel i: 

(c,g,i?)(i) = argmin ^ w{i - j) - (3) 

j—iG [—m.m] X [—rn.rn] 

We use different weight functions w with their support limited to the box size such 
as triangle, cosine or cosine^ tapers. 

The second method commonly used is to calculate the Taylor coefficients (1) not 
from the original but from a filtered image 

i{x) = j2Mx-j)nj) (4) 

i 

As window function Wd we use a normalised Gaussian of width d. The Taylor coef- 
ficients can now be explicitly derived by differentiation of 7, which however acts on 
the window function Wd instead on the image data. We therefore effectively use a 
filter kernel for each Taylor coefficient which relates to the respective derivatives of 
the window function Wd- 

One advantage of the latter method over the local fit described above is that the 
window width d can be chosen from M+ while the window size for the fit procedure 
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must be an odd integer 2m + 1. Both of the above methods regularize the Taylor 
coefficient estimate by the finite size of their window function. In fact, the window 
size could be considered as a regularization parameter. 

A common problem of regularized inversions is the proper choice of the regular- 
ization parameter. Lindeberg (1998) has devised a scheme for how this parameter 
can be optimally chosen. Our third method is a slightly modified implementation of 
his automated scale selection procedure. The idea is to apply method 2 above for 
each pixel repeatedly with increasing scales d and thereby obtain an approximation 
of the ridge's second derivative eigenvalues and , each as a function of the scale 
d. 

Since the h± and /ly are the principal second order derivatives of the image after 
being filtered with Wd, they depend on the width d of the filter window roughly in the 
following way. As long as d is much smaller than the intrinsic width of the ridge, d± 
= |(u"5^V)^log /|~^/^, the value in h± will be a (noisy) estimate of the true principal 
second derivative (tt^V)^/ of the image, independent of d. Hence, h± oc — Jmax/rfj_ 
for d^ <^ d]_. To reduce the noise and enhance the significance of the estimate, we 
would however like to choose d as large as possible. For d^ S> the result obtained 
for h± will reflect the shape of the window rather that of the width of the ridge. 

—Ima,xd±/d^ for d^ ':$> d\. Roughly, details near d ^ d± depend 
on the exact shape of the ridge, we have 

h±{d) =^^72 (5) 

For each pixel we consider in addition a quality function 

g(d) = dT(|/i^|-|/ijll), 7e(0,3) (6) 

which will vary as d'^ for small d and decrease asymptotically to zero for d ^ d± 
as d?~^. In between, q{d) will reach a maximum approximately where the window 
width d matches the local width d± of the ridge (which is smaller than the scale 
along the ridge). The choice of the right width d has now been replaced by a choice 
for the exponent 7. The result, however, is much less sensitive to variations in 7 
than to variations in d. Smaller values of 7 shift the maximum of q only slightly to 
smaller values of d and hence tend to favour more narrow loop structures. While 7 is 
a constant for the whole image in this automated scale selection scheme, the window 
width d is chosen individually for every pixel from the respective maximum of the 
quality factor q. 

In Figure 1 we show as an example a A = 171 A image of active region NOAA 
10930 observed by STEREO/SECCHI on 12 December 2006 at 20:43 UT and the 
corresponding image of q obtained with 7 = 0.75 and window sizes d in the range of 
0.6 to 4 pixels. Clearly, the q factor hass a maximiim in the vicinity of the loops. The 
distribution of the scales d for which the maximum q was found for each pixel is shown 
in Figure 2. About 1/3 of the pixels had optimal widths < one pixel, many of which 
originate from local longated noise and moss features of the image. The EUV moss is 
an amorphous emission which originates in the upper transition region (Berger et ai, 
1999) and are not associated with loops. For proper loop structures the optimum 
width found was about 1.5 pixels with, however, a widely spread distribution. 

In the case that i is located exactly on a ridge, u± is the direction across and it|| 
the direction along the ridge and h± and hn are the associated second derivatives 
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Figure 1. Original image of active region NOAA 10930 observed by STEREO/SECCHI 
(left) and the corresponding image of the quality factor q (6) obained from the third ridgel 
determination method by automated scale selection (right). This image was taken at A = 171 
A on 12 December 2006 at 20;43 UT and was not processed by the SECCHI prep routine. 
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Figure 2. Distribution of widths d of the window function for which (6) was found to 
reach a maximum when applied to the data in image (1). The maximum was independently 
determined for each image pixel for which h± < |. 



of the image intensity in the respective direction. A positive ridge is identified from 
the Taylor coefficients by means of the fofiowing conditions (Lindeberg, 1998) 

ttj^V/ = u^g = a vanishing gradient across the ridge (7) 
{u±V)^I = h± < a negative second order derivative (8) 



across the ridge 
a second order d 
or \h±\ > |/i||| across the ridge larger than along 



T 2 T 2 

l(itj^V) J| > |(it|| V) /| a second order derivative magnitude (9) 



The latter two inequalities are assumed to also hold in the near neighbourhood of 
the ridge and are used to indicate whether the pixel centre is close to a ridge. 
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Figure 3. Comparison of the resulting ridgcls from interpolation method (12, left) and (13, 
right). The images show an enlarged portion of tlic original data in Figure 1. The short sticks 
denote the local orientation of u± on the loop trace.) 



In the vicinity of the ridge, along a line x = i + u±t, the image intensity (1) then 
varies as 

I{t) ~ c + u±gt + u'±u±h±t^ (10) 

According to the first ridge criterion (7), the precise ridge position is where I{t) has 
its maximum. Hence the distance to the ridge is 

ulg 

(•max — ol, T ^ ' 

and a tangent section to the actual ridge curve closest to i is 

T 

r(s) = i "-^"^^ + suu for s G K (12) 

2h±u\^uj_ " 

Note that u'^u± = 1 for a unitary U. 

We have implemented two methods for the interpolation of the ridge position 
from the Taylor coefficients calculated at the pixel centres. One is the interpolation 

of the ridge centre with the help of (12). The second method interpolates the zeros 
of u^g in between neighbouring pixel centres i and j if its sign changes. Hence the 
alternative realization of (7) is 

if \c\ = \u±{i)u±{j)\ > Cmin 

and sign(c)(«Ig)(i) iu±g){i) < then 
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r = i + t{j - i) 



(13) 



where t = 



{g'^u±){i) - sign{c){gTuj_){j) 



The first condition insures that u±{i) and u±(j) are sufficiently parallel or an- 
tiparallel. Note that the orientation of uj_ of neighbouring pixels may be parallel or 
antiparallel because an eigenvector u± has no uniquely defined sign. 

In general, the interpolation according to (12) yields fewer ridge points along a 
loop but they have a fairly constant relative distance. With the second method (13), 
the ridge points can only be found at the intersections of the ridge with the grid lines 
connecting the pixel centres. For ridges directed obliquely to the grid, the distances 
between neighbouring ridge points produced may vary by some amount. Another 
disadvantage of the second method is that it cannot detect faint ridges which are 
just about one pixel wide. It needs at least two detected neighbouring pixels in the 
direction across the ridge to properly interpolate the precise ridge position. The 
advantage of the second method is that it does not make use of the second order 
derivative h± which unavoidably is more noisy than the first order derivative g. In 
Figure 3 we compare the ridgels obtained wth the two interpolation methods for the 
same image. 

The final implementation of identifying ridge points in the image comprises two 
steps: First the Taylor coefficients (1) are determined for every pixel and saved for 
those pixels wiiich have an intensity above a ttnesliold value /mim for which the 
ridge shape factor {h]_ — h^\)/ +h^\) exceeds a threshold Smin in accordance with 
(10) and which also satisfy (8). The second step is then to interpolate the precise 
sub-pixel ridge point position from tiie derivatives at these pixel centres by either of 
the above methodes. This interpolation complies with the third ridge criterion (7). 
The ridgel orientation u± are also interpolated from the cell centres to the ridgel 
position. The information retained for every ridge point n in the end consists of its 
location r„, and the ridge normal orientation tt_L,n defined modulo tt. 

2.2. Ridgel Connection to Chains 

The connection algorithm we apply to the ridgels to form chains of ridgels is based 
on the cocircularity condition of Parent and Zucker (1989). 

For two ridgels at r„ and r„+i a virtual centre of curvature can be defined which 
forms an isoceles triangle with the ridgels as shown in Figure 2.2. One edge is formed 
by the connection between the two ridgels of mutual distance hnM+i - The two other 
triangle edges in this construction connect one of the two ridgels with the centre of 
curvature which is chosen so that these two symmetric edges of the isoceles triangle 
make angles Aa„,„+i and Aan+i,n as small as possible with the respective ridgel 
orientation it_L,„ and it_L,n+i, respectively. It can be shown that 



requires equal magnitudes for the angles Aan,n+i and Aan+i,n- The distance rn,n+i 
= i'n+i,n is the local radius of curvature and can be calculated from 



min (Aa„,„+i + Aa„+i_„) 



(14) 




(15) 



cos {^{an,n+l + (Xn+l,n)) 
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Figure 4. Illustration of angles and distances of a connection element between two ridgels 
according to the cocircularity condition of Parent and Zucker (1989). The ridgel positions are 
indicated by a small dot, the ridge normal orientation by the line centred at the ridgel position. 
The axis units are image pixels. 

where an,n+i is the angle between rn+i — Vn and ±u_l_„, the sign being chosen so 
that |aTi,n+i| < 7r/2. 

With each connection between a pair of ridgels we associate a binding energy 
which depends on the parameters derived above in the form: 

e„,n+i = (— ) +[- ) +(-r- ) -3 (16) 

(-Kmax ' n,n+l 'imax 

Note that Aa^ ^+i = ^cen+i,n according to the cocircularity construction and hence 
£n,n+i is symmetric in its indices. The three terms measure three different types of 
distortions and can be looked upon as the energy of an elastic line element. The 
first term measures the deviation of the ridgel orientation from strict cocircularity, 
the second the bending of the line element, and the third term its stretching. The 
constants amax, "rmim and ft-max give us control on the relative weight of the three 
terms. rmi„ is the smallest acceptable radius of curvature and /imax is the largest 
acceptable distance if we only accept connections with a negative value for the energy 
(16). 

In practical applications, the energy (16) is problematic since it puts nearby ridgel 
pairs with small distances hn,n+i at a severe disadvantage because small changes of 
their u± easily reduces the radius of curvature r„,„+i below acceptable values. We 
therefore allow for measurement errors in r and u± and the final energy considered 
is the minimum of (16) within these given error bounds. 

The final goal is to establish a whole set of connections between as many ridgels 
as possible so that the the individual connections add up to chains. Note that each 
ridgel has two "sides" defined by the two half spaces which are separated by the 
ridgel orientation u^. We only allow at most one connection per ridgel in each of 
these "sides" . This restriction avoids junctions in the chains that we are going to 
generate. The sum of the binding energies (16) of all accepted connections ideally 
should attain a global minimum in the sense that any alternative set of connections 
which complies with the above restriction should yield a larger energy sum. 
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Figure 5. Sketch of the curve-fit parameters. The ridgels are represented by their location and 
the two-pixel-long bar of the rpc orientation. For each ridgel i, the proximity to the smooth-fit 
curve is expressed by their distance and the angle Act; between the distance direction of di 
to the curve and the rpc orientation. Another measure of the quality of the curve is the inverse 
radius of curvature r. 

We use the following approach to find a state which comes close to this global 
minimum. The energy en,n+i is calculated for each ridgel pair less than /imax apart, 
and those connections which have a negative binding energy are stored. These latter 
are the only connections which we expect to contribute to the energy minimum. Next 
we order the stored connections according to their energy and connect the ridgels 
to chains starting from the lowest-energy connection. Connections to one side of 
a ridgel which has already been occupied by a lower energy connection before are 
simply discarded. 

2.3. Curve Fits to the Ridgel Chains 

In this final section we calculate a smooth fit to the chains of ridgels obtained above. 
The fit curve should level out small errors in the position and orientation of individual 
ridgels. We found from experiments that higher-order spline functions are far too 
flexible for the curves we aim at. We expect that magnetic-field lines in the corona 
do not rapidly vary their curvature along their length and we assume this also holds 
for their projections on EUV images. We found that parametric polynomials of third 
or fifth degree are sufficient for our purposes. Hence for each chain of ridgels we seek 
polynomial coefficients Qn which generate a two-dimensional curve 

5 

p(i) = ^g„r for iG[-l,l] (17) 

which best approximates the chain of ridgels. What we mean by "best approx- 
imation" will be defined more precisely below. The relevant parameters of this 
approximation are sketched in Figure 2.3. 
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Figure 6. Fit curves obtained for those chains which involve ten or more ridgels (left) and 
curves remaining after cleaning of those curves which are due to moss (right). 



The polynomial coefficients q of a fit (17) are determined by minimising 
J2 {dfd^)+fi{p"'^p"){ti) 

zGchaiii 

wliere di = ri - p{ti) 



(18) 



with respect to q„ for a given /i. Initially, we distribute the curve parameters ti in the 
interval [—1, 1] such that the differences ti — tj of neighbouring ridgels are propor- 
tional to the geometric distances \ri — rj\. The p" are the second-order derivatives 
of (17). Hence, the second term increases with increasing curvature of the fit while 
a more strongly curved fit is required to reduce the distances di between and the 
first order closest curve point p(ti). 

The minimum coefficients qnifJ') can be found analytically in a straight forward 
way. Whenever a new set of qnifJ-) has been calculated, the curve nodes U are 
readjusted by 



argmin((r,; — p{t)y 



(19) 



so that p{ti) is always the point along the curve closest to the ridgel. 

For different fi this minimum yields fit curves with different levels of curvature. 
The local inverse radius of curvature can at any point along the curve be calculated 
from (17) by 



1 _ \p"{t)xp\t)\ 
r{t) |p'(t)|3 

The final n is then chosen so that 



(20) 



Echa 



.(M)= E 

i^chain 



f^max 



+ E 

I G chain 



Q^max 



+ r„ 



I r{tf 



dt 



(21) 
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is a minimum where Aan is the angle between the local normal direction of dn of 
the fit curve and the ridgel orientation ±ttj^_„, the sign again chosen to yield the 
smallest possible |Aaj| The meaning of the terms is obvious, and clearly the first 
two terms in general require a large curvature which is limited by the minimization 

of the last term. 

Expression (21) depends nonlinearly on the parameter fx which we use to control 
the overall curvature. The minmum for (21) is found by iterating /n starting from 
a large numerical value, i. e., Qb straight line fit. Xhe parameters /"mini cKmax? and 
dmax can be used to obtain fits with a different balance between the mean square 
spatial and angular deviation of the fit form the "observed" chain of ridgels and the 
curvature of the fit. Unless these parameters are chosen reasonably, e.g. rmin not too 
small, we have always found a minimum for (21) after a few iteration steps. 

In the left part of Figure 6 we show the final fits obtained. For this result, the 
ridgels were found by automated scaling and interpolated by method (12), the param- 
eters in (16) and (21) were hmax = 3.0 pixels, rmin = 15.0 pixels, and amax = 10.0 
degrees. The fits are represented by fifth degree parametric polynomials. 

Obviously, the image processing cannot easily distinguish between structures 
which are due to moss and bright surface features and coronal loops. Even the 
observer is sometimes misled and there are no rigorous criteria for this distinction. 
Roughly, coronal loops produce longer and smoother fit curves, but there is no strict 
threshold because it may appear that the fit curve is split along a loop where the 
loop signal becomes faint. As a rule of thumb, a restriction to smaller curvature by 
choosing a higher parameter rmax and discarding shorter fit curves tends to favour 
coronal loops. Eventually, however, loops are suppressed, too. We have therefore 
appended a user interactive tool as the last step of our processing which allows us 
to eliminate unwanted curves, merge or split curves when smooth fits result with an 
energy (21) of the output fits not much higher than the energy of the input. The left 
part of Figure 6 shows the result of such a cleaning step. 

3. Application 

In this section we present an application of our segmentation tool to another EUV 
image of active region NOAA 10930 taken by the SECCHI instrument onboard 
STEREO A. This EUV image was observed at A = 195 A on 12 December 2006 at 
23:43:11 UT. At that time the STEREO spacecraft were still close so that stereoscopy 
could not be applied. We therefore selected an image which was taken close to 
the MDI magnetogram observed at 23:43:30 UT on the same day. It is therefore 
possible to calculate magnetic-field lines from an extrapolation model and project 
them onto the STEREO view direction to compare them with the loop fits obtained 
with our tool. In Figure 7 the MDI contour lines of the line-of-sight field intensity 
were superposed on the EUV image. 

The loop fits here were obtained by applying the automated scaling with d up to 
two pixels, i.e. window sizes up to 2d -f 1 = 5 pixels, to identify the ridgels. Pixels 
with maximum quality q below 0.4 were discarded and we applied method (12) to 
interpolate the local ridge maxima, ftmax = 5.0 pixels, rmin = 15.0 pixels, and 
flmax = 10.0 degrees. The fits are fifth degree parametric polynomials. In Figure 8, 
we show some of the fits obtained which are most likely associated with a coronal 
loop. They are superposed onto the EUV image as red lines. Those loops which were 
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Figure 7. MDI contours overlaid on the STEREO/SECCHI EUV image for NOAA 10930. 
The EUV and MDI data were recorded on 12 December 2006 at 23:43:11 UT and 23:43:30 
UT, respectively. The colour code on the right indicates the field strength at the contour lines 
in Gauss. 

found close to computed magnetic-field lines are displayed again in the left part of 
Figure 9 with loop numbers so that they can be identified. 

The magnetic-field lines were computed from the MDI data by an extrapolation 
based on a linear force-free field model (see Seehafer, 1978 and Alissandrakis, 1981 for 
details). This model is a simplification of the general nonlinear force-free magnetic- 
field model 

V X B = qB where B • Va = (22) 

and a may vary on different field lines. An extrapolation of magnetic surface obser- 
vations based on this model requires boundary data from a vector magnetograph. 
The linear force-free field model treats a as a global constant. The advantage of 
linear the force-free field model is that it requires only a line-of-sight magnetogram, 
such as MDI data, as input. 

A test of the validity of the linear forec-free assumption is to determine different 
values of a from a comparison of field lines with individual observed loop structures 
(e.^., Carcedo et al., , 2003). The range of a obtained then indicates how close the 
magnetic field can be described by the linear model. Since the linear force-free field 
has the minimum energy for given normal magnetic boundary field and magnetic 
helicity, a linear force-free field is supposed to be much more stable than the more 
general non-linear field configuration (Taylor, 1974). 

We calculated about 5000 field lines with the linear force-free model with the a 
value varied in the range from -0.0427 Mm~^ to 0.0356 Mm~^ These field lines were 
then projected onto the EUV image for a comparison with the detected loops. 

For each coronal loop li, we calculate the average distance of the loop to every 
projected field line hj discarding those field lines that do not fully cover the observed 
loop li. This distance is denoted by Ci.{bj). For details of this distance calculation see 
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Figure 8. The left diagram shows as red lines the loops identified by the segmentation tool 
for the EUV image in Figure 7. In order to show the loops more clearly, the image in the 
background was contrast enhanced by an unsharp mask filter. The right diagram displays in 
yellow field lines calculated from the MDI data. The fieldlines were selected so that they are 
located closest to the extracted loops in the left part of the image. See the text for more details 
on how the field lines were determined. 
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Figure 9. The left panel shows the loops identified by the segmentation tool which corre- 
sponds to the closed magnetic filed lines; In the right panel are the loops (solid lines) with 
their best-fitting field lines(dotted lines), x and y axis are in units of EUV pixels. 



Feng et al, (2007). In the end we could find a closest field line for every coronal loop 
by minimizing Ci-(bj). The detected loops and their closest field lines are plotted in 
the right diagram of Figure 9. An overplot of the closest field lines onto the EUV 
image is shown in Figure 8 

In Table 1 we list the distance measure C along with the loop number and the 
linear force-free parameter a for the closest field line found. We find that our a 
values are not uniform over this active region, that is, the linear force-free model is 
not adequate to describe the magnetic properties of this active region. This is also 
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Table 1. Identified loops, a values of the best 
fitting field lines and the averaged distances in 
units of pixel between the loop and the best 
fitting field line. 



Loop No. 


C((b) (pixel) 


a(10-3Mm-i) 


1 


3.1957 


-10.680 


3 


3.7523 


-35.600 


4 


2.3432 


-9.2560 


11 


10.7692 


-35.600 


13 


0.4864 


2.1360 


14 


1.3636 


-9.2560 


15 


4.2386 


17.088 


17 


4.8912 


16.376 


18 


2.4256 


-14.240 


19 


2.5388 


-32.752 



seen by the characteristic deviation at their the upper right end in Figure 8 between 

the eastwards-inclined loops (solid) and their closest, projected, field lines (dotted). 
With no value of a the shape of these loops could be satisfactorily fitted. Further 
evidence for strong and inhoinogeneous currents in the active region loops is provided 
by the fact that only 2.5 hours later, at 02:14 UT on 13 December, a flare occurred 
in this active region and involved the magnetic structures associated with loops 2, 
4, and 17. 



4. Discussion 

EUV images display a wealth of structures and there is an important need to reduce 
this information for specified analyses. For the study of the coronal magnetic field, 
the extraction of loops from EUV images is a particularly important task. Our tool 
intends to improve earlier work in this direction. Whether we have achieved this 
goal can only be decided from a rigourous comparison which is underway elsewhere 
(Aschwanden et al., 2007). At least from a methodological point of view, we expect 
that our tool should yield improved results compared to Strous (2002, unpublished) 
and Lee, Newman, and Gary (2006). 

From the EUV image alone it is often difficult to decide which of the features 
are associated with coronal loops and which are due to moss or other bright surface 
structures. A final comparison of the loops with the extrapolated magnetic field 
and its field line shapes is therefore very helpful for this distinction. Yet we have 
avoided to involve the magnetic-field information in the segmentation procedure 
which extracts the the loops from the EUV image because this might bias the loop 
shapes obtained. 

For the case we have investigated, we find a notable variation of the optimal a 
values and also characteristic deviation of the loop shapes from the calculated field 
lines. These differences are evidence of the fact that the true coronal magnetic field 
near this active region is not close to a linear force-free state. This is in agreement 
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with earlier findings. Wiegelmann et ai, (2005), e.g., have shown for another active 
region that a nonlinear force-free model describes the coronal magnetic field more 
accurately than linear models. The computation of nonlinear models is, however, 
more involved due to the nonlinearity of the mathematical equations {e.g., Wiegel- 
mann, 2004; Inhester and Wiegelmann, 2006). Furthermore, these models require 
photospheric vector magnetograms as input, which were not available for the active 
region investigated. 

Coronal loops systems are often very complex. In order to acess them in 3D, the 
new STEREO/SECCHI telescopes now provides EUV images which can be analysed 
with stereoscopic tools. We plan to apply our loop-extraction program to EUV images 
from different viewpoints and undertake a stereoscopic reconstruction of the true 3D 
structure of coronal loops along the lines described by Inhester (2006) and Feng 
et ai, (2007). The knowledge of the 3D geometry of a loop allows to estimate more 
precisely its local EUV emissivity. From this quantity we hope to be able to derive 
more reliably the plasma parameters along the length of the loop. 

Other applications can be envisaged. An interesting application of our tool, e.g., 
will be the investigation of loop oscillations. Here, the segmentation tool will be 
applied to times series of EUV images. We are confident that oscillation modes and 
in the case of a STEREO/SECCHI pairwise image sequence the polarisation of the 
loop oscillation can also be discerned. 
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